Semiclassical Wigner distribution for two-mode entangled state 
generated by an optical parametric oscillator 



0^ 
O 
O 
(N 



K. DechoumE M. D. Hahn, and A. Z. KhourjH 

Instituto de Fisica, Universidade Federal Fluminense, 
Av. Gal. Milton Tavares de Souza s/n, 24210-346 Niteroi - RJ, Brazil 

R. O. Vallejoil!] 

Centra Brasileiro de Pesquisas Fisicas, 

Rua Dr. Xavier Sigaud 150, 
22290-180 Rio de Janeiro - RJ, Brazil 
(Dated: June 29, 2009) 

We derive the steady state solution of the Fokker-Planck equation that describes the dynamics 
of the nondegenerate optical parametric oscillator in the truncated Wigner representation of the 
density operator. The adiabatic limit of strong pump damping is assumed. This phase space image 
provides a clear view of the intracavity two-mode entangled state valid in all operating regimes of 
the OPO. A nongaussian distribution is obtained for the above threshold solution. 
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I. INTRODUCTION 

Among the quasiprobability functions that represent 
the density operator of a quantum state, the Wigner dis- 
tribution has undoubtedly advantages over others, since 
in this phase space representation the quantum-classical 
correspondence is, in general, much more visible. At the 
same time these functions contain all information avail- 
able in the density operator. 

Unfortunately, just a few states have been repre- 
sented by this distribution due to the difficulty to solve 
Fokker-Planck equations for nonlinear systems. Indeed, 
the phase space description of quantum systems is well 
known for quadratic hamiltonians, but very little is 
known for nonlinear systems, such as, for example, the 
single-mode degenerate parametric oscillator Q, or the 
transverse multimode degenerate parametric oscillator 

i- 

Of special interest is the two-mode entangled state gen- 
erated in the optical parametric oscillator (OPO), used in 
many experiments of quantum information. The Wigner 
distribution of this state was obtained exactly using the 
solution derived from the complex P-representation Q, 
and tranformed into a Wigner distribution. Since this 
analytical result appears as an infinite series of gamma 
functions (hypergeometric series), some important phys- 
ical aspects are hidden by its mathematical complexity. 
For example, two-mode entanglement is not of easy iden- 
tification from the usual partial transpose criterion [J Q . 
The same happens to the generalized P-representation for 
the intracavity parametric oscillator derived in Ref.j^. 
In order to calculate the intensity correlations, the au- 
thors obtained an integral expression that depends on 
degenerate hypergeometric functions, and the outcome 
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also depends on an infinite sum. 

In the present work, we derive the Wigner function of 
a two-mode quantum state of the electromagnetic field 
generated by an OPO using a truncated Wigner equa- 
tion. This truncated approximation is valid when the 
nonlinear! ty is small enough, which is true for most of 
the OPO experiments developed so far. A simple expres- 
sion is obtained for the Wigner distribution, which ren- 
ders clear the many quantum features of the intracavity 
OPO field such as squeezing and entanglement. More- 
over, this expression is valid in all regions of the OPO 
operation regimes, that is, below, and above threshold, 
where the validity of both the linear approximation, and 
the perturbation theory break down. This distribution 
also describes the OPO quantum fluctuations properly 
around the operation threshold, while the usual linear 
approach provides divergent results. In order to illus- 
trate the reliability of the Wigner function we obtain, 
moments calculated from this distribution are compared 
with the exact quantum results, for any values of the 
control parameters of the OPO. 



II. THEORETICAL MODEL 

We present here a model of three quantized modes cou- 
pled by a nonlinear crystal inside a triply resonant Fabry- 
Perot cavity. The Heisenberg picture Hamiltonian that 
describes this open system is given byi?]] 



H 



+ ih(Ee 



a a. 



J2 («^rl + a]r. 



1=0 



(1) 



Here E represents the external coherent driving pump 
field at frequency loq. The operators oq, oi and 0,2 rep- 
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resent the pump, signal and idler fields, respectively, 
satisfying the following frequency matching condition, 
luq — LOi + UJ2. The terms Ti represent damping reser- 
voir operators, and x is the nonlinear coupling constant 
due to the second order polarizability of the nonlinear 
crystal. 

The master equation for the reduced density operator, 
after the elimination of the heat bath by standard tech- 
niques ^7] , is given by 
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where 7^ is the corresponding mode damping rate. 

In order to treat the operators evolution, we now turn 
to the method of operator representation theory. These 
techniques can be used to transform the density ma- 
trix equation of motion into c-number Fokker-Planck or 
stochastic equations. We can write down the master 
equation in the Wigner representation, by using the fol- 
lowing characteristic function 
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so that the Wigner distribution can be written as Fourier 
transform of the characteristic function: 
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The phase space Wigner equation for the nondegener- 
ate parametric amplifier that corresponds to the master 
equation ^ is then 
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This is not a Fokker-Planck equation due to the third 
order derivative term, but in the case where x is small 
enough we can drop this term. The truncated equation 
so obtained is a genuine Fokker-Planck equation with a 
positive diffusion term, and we can easily derive the cor- 
responding stochastic differential equations in the rotat- 
ing frame {aj — aj exp {—iujjt) with j = 0, 1, 2): 
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It is possible to find a stationary solution of the truncated 
Fokker-Planck equation when we adiabatically eliminate 
the pump mode variables. This means that we are con- 
sidering the relaxation rate of this mode much larger than 
those of the downconverted modes, that is 70 3> 71,72- 
The stationary solution for the pumped mode is 

ao = — [-B - X'5i<52 + ^/7oCQ(^)] , (7) 
70 

where the noise term is retained in the adiabatic elimi- 
nation in order to properly deal with the noise dynamics. 
We are then left with two nonlinear dynamical equations 
for the complex amplitudes of the down-converted fields, 
where the pump amplitude is replaced by expression ([7]). 
We now define the following real quadrature variables: 
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so that the remaining dynamical equations can be cast 
in the compact form: 
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where X = [xi,yi,X2,y2]'^ is a column vector, and the 
drift matrix is defined as 
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We have set 71 =72 =7, which is a reasonable physical 
assumption for most OPO experiments, and defined the 
normalized pump parameter p = X^/(77o)j as well as 
the nonlinear coupling g = xl (\/277o)- B is a 4 x 6 
matrix defined as 
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and ^(t) is a six component column vector whose entries 
are uncorrelated real Gaussian noises associated with the 
noise terms for the three interacting fields. 

Now this set of stochastic equations can be mapped 
onto a genuine Fokker-Planck equation: 
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where the diffusion matrix is defined as D = BB"^, and 
summation over repeated indices is assumed throughout 
the text. 

This Fokker-Planck equation admits a steady state so- 
lution in the potential form VF(X) = A^exp (— / dXi) 
with Zi given by 



2A, 



(13) 



This solution provides a quite simple form for the 
Wigner distribution which allows for the calculation of 
several statistical properties. It is important to no- 
tice that the potential solution is achievable only when 
7i = 72 = 7 as we have already assumed [l^ . In this 
case, the steady state Wigner distribution for signal and 
idler reads: 
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where J\f is the normalization constant and unimportant 
terms 0{g'^) were neglected in the argument of the po- 
tential solution. 

As a concrete physical example, expression (I14p can 
describe the joint Wigner distribution for two polariza- 
tion modes of a frequency degenerate type II OPO, where 
signal and idler are distinguished by their polarization 
states. All statistical properties, including experimen- 
tally accessible quantities like quadrature noise and cor- 
relations in the stationary state may be calculated with 
this distribution in any OPO operation regime. It is 
interesting to notice that the distribution given by ex- 
pression is single peaked below threshold operation 
(/i < 1) and becomes double peaked above threshold. 
In Fig. (1), it is possibible to visualize the conditional 
Wigner distribution for fixed values of yi and 2/2, for an 
OPO operating below, at, and above threshold. It is easy 
to identify the squeezed and unsqueezed quadratures in 
all regimes. 

For the nonlinear media employed in most OPO de- 
vices, the coupling parameter g <^ 1. Therefore, the 
first two terms in the exponent of the Wigner distribu- 
tion governs the OPO behaviour below threshold. They 
are enough to predict the usual nonclassical features. In 
fact, it is evident now that the usual linearized approach 
to quantum noise in below-threshold OPOs is equivalent 





Figure 1: Conditional Wigner distribution for j/i = and 
y2 — 0, a) below (/i = 0.5); b) at (/i = 1.0); and c) above 
threshold = 1.5). In all cases = 0.01. 



to completely neglect the g^ term in the distribution. 
However, it is important to notice that for /i > 1, the 
Wigner distribution becomes divergent if we neglect the 
g^ term. This means that this term plays a crucial role 
for a complete description of the OPO behaviour at and 
above threshold. 

Additional insight is provided by a careful investiga- 
tion of the marginal distribution obtained by integrating 
W{xi,yi,X2,y2) with respect to one mode variables. Of 
course, due to the symmetry of W, the form of the re- 
sulting marginal distribution is independent of the mode 
being traced out. The marginal distribution for mode 2 
becomes: 



W{x2,y2) = 27rA/'x 



(15) 
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It gives the statistical properties of isolated measure- 
ments on mode 2. The variances given by this marginal 
distribution are larger than those of the vacuum state. 
This excess noise can be seen as a consequence of informa- 
tion loss when one looks only at part of the whole system. 
Fig. (2) presents the marginal distribution below, at and 
above threshold. Above threshold the marginal distri- 
bution is clearly peaked out of the phase space origin, 
what is a consequence of the macroscopic amplification 
of the down-converted fields. However, the distribution 
presents radial symmetry indicating complete phase un- 
certainty in each individual mode. 






Figure 2: Marginal distributions: a) below threshold, ^ = 0.8; 
b) at threshold, fj, — 1.0; c) above threshold, p = 1.2. In all 



cases g 



0.01. 



III. COMPARISON WITH THE LINEARIZED 
QUANTUM APPROACH 

Almost the totality of theoretical works devoted to the 
analysis of quantum noise in optical systems rely on lin- 
earization of the small quantum fluctuations around the 
macroscopic steady state mean values. It is worth to no- 
tice that this procedure has limited validity, especially 
around the threshold critical point where quantum fluc- 
tuations may become comparable to the mean values. In 
order to evidence the breakdown of the linearized theory, 
we shall compare its results with those obtained from the 
solution of the Fokker-Planck equation (fT2|) . 

The usual approach given to quantum noise in the 
literature is obtained from first order stochastic equa- 
tions of motion [l^, ■ These equations are often used 
to predict squeezing in a linearized fluctuation analy- 
sis. They are non-classical in the sense that they can 
describe states without a positive Glauber-Sudarshan P- 
distribution, but correspond to a Gaussian Wigner dis- 
tribution. 

We now find it useful to introduce combined field 
quadratures, as in two- mode approaches used previ- 
ously [To| . These combined quadratures are the Einstein- 
Podolsky-Rosen (EPR) variables used to characterize 
continuous variable entanglement between the down con- 
verted fields. They are defined as 
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These quantities correspond to the squeezed and anti- 
squeezed combined quadratures obtained in the lin- 
earized theory. /^From the first order stochastic equa- 
tions one can easily obtain the steady state variances of 
the EPR variables [11, 111: 
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for below threshold operation, and 
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for above threshold operation. 

Now let us briefly discuss the predictions of the lin- 
earized approach and its validity. The quadratures 
and y+ exhibit the expected squeezing, going from the 
vacuum fluctuations for zero pump until 50% intracavity 
squeezing (which corresponds to perfect squeezing out- 
side the cavity) at and above the oscillation threshold. 
The unsqueezed quadratures x^ and y_ present diver- 
gent behaviour around threshold, which is certainly not 
physical. In fact, as we shall see next, the steady state 
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solution of the Fokker-Planck equation (fT2|) gives a well 
behaved Wigner distribution which does not display any 
divergences. This Wigner distribution will also put limits 
on the amount of squeezing attainable. 

In order to investigate the two-mode entanglement di- 
rectly from the Wigner distribution, we now write it in 
terms of the EPR variables. Below threshold, where the 
linearization of the equation of motion is valid, we can 
neglect the term, and the distribution can be approx- 
imated as 

WL{x+,y+,x_,y_) =7Vea;p|-i [(1 + ^) x^_ + 

(1 + m) y\ + {l- /i) x\ + {l- m) yl] } . (19) 

With this expression we can easily calculate any mo- 
ment of the distribution. For instance, we easily reobtain 
the results of Eqs. p7p for the intracavity noise squeezing 
in the combined quadratures. As we mentioned before, 
at threshold we have perfect external squeezing while the 
unsqueezed combined quadratures blow up showing the 
failure of linearization [l^ . Moreover, it is easy to char- 
acterize this state as entangled by using the Duan-Simon 
criterion [1, . Above threshold we can also calculate the 
moments and we find that this criterion shows the sufi- 
cient condition to characterize this state as entangled, in 
spite of not being a Gaussian state in that regime. 

It is also interesting to compare the results obtained for 
{x\) with the linearized theory (both below and above 
threshold) with those obtained from the Wigner func- 
tion given by Eq.lHl]). This is done in Fig.®. While 
the linearized approach gives a divergent behavior at the 
oscillation threshold, the results obtained from the trun- 
cated Wigner function gives a more realistic smooth be- 
havior in agreement with the full quantum solution of 
Refs.jl, @, [l^. The results for the squeezed quadrature 
(x^) are shown in Fig|3]where, again, the results given by 
the truncated Wigner distribution agree with those ob- 
tained from the full quantum solution. Below threshold, 
the linearized approach also agrees with the full quan- 
tum theory and the truncated Wigner approach, but 
an important deviation can be observed above thresh- 
old where the linearized theory gives (a;^) = 0.5 for any 
pump power. The same results are obtained for the other 
squeezed variable so that the Duan-Simon separa- 

bility criterion 0, Q shows that the two-mode quantum 
state is entangled. However, it is important to notice 
that both the full quantum approach and the truncated 
Wigner distribution predict a less pronounced violation 
of the criterion above threshold. 

In Refs.jl, [l^ Fokker-Planck equations were derived 
for the complex-P ^ and positive-P [l5| representations, 
and stationary solutions were obtained for the corre- 
sponding distributions. Since these Fokker-Planck equa- 
tions were derived without any truncation, we can con- 
sider the stationary distributions so obtained as exact. In 
Ref. 3] an exact stationary Wigner function was mapped 
from the complex-P distribution of Ref. and expressed 
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Figure 3: Comparison between the values of the anti-squeezed 
variable {x\) obtained from the linearized theory (dashed 
line) with those obtained from the Wigner function of ref- 
eren Q ' ^ ' ■ — 1" '[ ' ;s) . 

In a 





1 1 1 1 








u + + ^ * * * 


1 . . . 1 . . . 



Figure 4: Comparison between the values of the squeezed 
variable (a;!) obtained from the linearized theory (solid line) 
with those obtained from the Wigner function of Ref. [3] 
(circles), and the one given by Ea. (|14p (crosses). In all cases, 
we have set — 0.01. 



in terms of Bessel functions. It is important to remark 
that the simple Wigner function presented here provides 
a very good approximation for those exact solutions, as 
can be seen in Figsl3]and|4l even for a rather large value 
for the nonlinear coupling g^. 



IV. CONCLUSION 

By truncating the evolution equation for the two-mode 
Wigner function of an optical parametric oscillator, we 
derived a Fokker-Planck equation which admits a sim- 
ple potential solution in any operation regime of the 
OPO. For the unsqueezed quadratures, this solution does 
not present the unphysical divergent behavior of the lin- 
earized theory on the oscillation threshold. Moreover, 
it provides good agreement with exact quasiprobabilities 
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already available in the literature, both for the squeezed 
and unsqueezed quadratures. While these exact solu- 
tions are given as infinite series, the potential solution of 
the truncated Wigner distribution presents a simple form 
allowing for an easier visualization of the phase space dis- 
tribution. 
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